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VII. Direct and expeditious methods of calculating the excentric 
from the mean anomaly of a planet. By the Reverend Abram 
Robertson, D. D. F.R. S. Savilian Professor of Astronomy in 
the University of Oxford, and Radcliffian Observer. Commu- 
nicated by the Right Hon. Sir Joseph Banks, Bart. G. C. B. 

P. R. S. 

Read February 15, 1816. 

Oince the publication of Kepler's discoveries in astronomy, 
the attention of men of science has frequently been directed 
to the problem distinguished by his name, and their exertions 
have frequently been employed to overcome the acknowledged 
difficulty of its solution. A statement of the various degrees 
of success, with which these endeavours have been made, is 
foreign to the present design. An account of this kind is now 
also needless, as Dr. Brinkley's examination of such attempts, 
published in the ninth volume of the Transactions of the Royal 
Irish Academy, affords a satisfactory review of most of the 
proceedings on this subject, previous to the year 1802. 

After the following methods had occurred to my consider- 
ation, and I had fully proved their utility by actual application 
to examples, I was anxious to ascertain whether any author 
had anticipated me in the manner in which the investigations 
are conducted. With this view I examined such solutions as 
are referred to in Dr. Brinkley's very able Memoir, all those 
mentioned by Montucla,* of which I could procure a sight, 

# Histoire des Mathematiques, Tom. II. p. 343, &c. I have searched, without 
success, for Lorgna's and Trembley's publications. 
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and some others which had occurred to me in the course of 
my reading. The result of this search is a belief that no one 
before has aimed at a direct solution through the same small 
angle, and, by means of an equation in which this angle and its 
powers are the only unknown quantities, obtained a quickly 
converging series for its value in known terms. The small. 
angle being found, with due precision, the eccentric anomaly 
readily becomes known. 

M. Delambre, in his Astronomy, published at Paris in 
3814, in three quarto volumes, calculates the excentric ano- 
maly by a method founded on those of Cassini,* La Caille,-|" 
Simpson,^ and Cagnoli,§ This eminent astronomer says of 
it, " Ce procede, le plus clirecte que je connaises, est aussi le 
plus precis ; il n'est qu'approximatif, mais il est toujours 
exact au-dela des dixiemes de seconde pour toutes les pianetes 
de notre systeme." That the reader may readily judge how 
far the third method, which I now propose/ deserves atten- 
tion, I have annexed to my investigation the calculation by it 
of M. Delambre's two examples, and also one relating to the 
comet of 1682 and 1759. 

Each of the following methods of solution is to be considered 
as direct, although it proceeds through the medium of what 
is commonly called Cassini's approximation. This approxi- 
mation, as here used, can only be considered as the first cer- 
tain step in the computation. No hypothesis is introduced 

* Memoires de l'Academie de 1'annee 1719. 

f Le^ns Elementaires d'Astronomie, Paris, 1761. 

% Essays on several curious and useful subjects, London, 1740. 

§ Trigonometries Paris, 1786. 
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into the proceeding, and therefore no correction by trial and 
error is requisite. 




Let ALP be the orbit of a planet, C the centre of the ellipse, 
S that focus in which the sun is placed, and AMP a circle 
described on the greater axis AP as a diameter. Let L be 
the true place of the planet, and AM the corresponding mean 
anomaly. Through L let the straight line ER be drawn, 
perpendicular to AP, and let it meet AP in R and the circle in 
E. Let EC, CM, SM be drawn, and let ST be perpendicular 
to EC, SG to CM and MN to ST. Then MN is parallel to 
ET, and NT is equal to the sine of the arc EM. It is easily 
proved, as in almost every writer on the subject, that ST is 
equal to the arc EM. 

In this Problem it is supposed that AC, CS, AM are given, 
and it is required to find AE the excentric anomaly, for AE 
being found, the true anomaly ASL is easily obtained. 

In each of the three following methods the angles CMS, 
CSM are used, and their difference is found by this pro- 
portion. 

mdcccxvi. S 
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CM + CS : CM— CS : : tan . \ ACM : tan. \ (CSM — CMS). 
Hence the angles become known by their sum and difference. 
As the angle SMN is very small, and consequently the angle 
MSC nearly equal to EGA in the orbits of almost all the 
planets, this way of finding the angle CSM is usually called 
Cassini's approximation to the excentric anomaly ACE. 

First Method. 

Having found the angle CSM, sin. CSM : CM : : sin. SCM : 
SM, which therefore becomes known. Let % equal the angle 
SMN, s equal the series expressing its sine, and c equal the 
series expressing its cosine. Put a equal the sine of CMS, 
and h equal its cosine. Then, radius being i , ac — bs = sin. 
CMN = sin. ECM = sin. (CMS —z). Also, 1 : SM : : s : 
SN = SM x s, and ac-fa+SMx s = TN + SN = EM = 
CMS — %, and therefore CMS = % -f ac — bs + SM * s = % -\> ac 



+ SM — bxs. 
Letd=SM— 6,andthenCMS=:« + flc + A=« + a(i-r- ~ 

+ ^1 « + *. _&c.)+<*(*— -+ — &c!> 

"2.3.4 23.4.5.6 ' 2.3.4.5.6.7.8 / 1 V 2.3 • 2.3.4.5 ; 

+ , , a& z dz 2 . ax.* , d% ? 

~ 2 2.3 ' 2.34 ' 2.3.4.5 

Let «=CMS~a, and putting A, B, C, &c. for the coefficients 

e — Ax - B% 2 — . Cz 3 + D% 4 + Ez 5 — F% 6 — G% 7 + &c. 

By reversing this series we find % = — -— , e -j- -rrrrr^ e*' 

•' ° 1 + a ' 2(14-") 

L - 3 . a* 3 a 4 i 5^ 4 1 S fl3 4. 

"TeTT+d^ 6 "r 2 (i+d)5 e 24(1 +d) s "t~52(i+d)« e T 8[i+d)»* 

-f- &C. 

This equation is in parts of the radius, and in order to have 
it in degrees we use this proportion, 1 : 57 , 2957795 :.;z: S^ 
2957795 * = R ° %> putting R° for 57 . 2957795. 
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Hence R° z = ^ e + ^™ e 2 + ^-, e 3 + —^ e 3 + &c. 

Second Method. 

The substitutions for the sines and cosines of the angles 
CMS, SMN being as in the preceding method, let SG be per- 
pendicular to MC, and then radius being 1, 1 : CS : : sin. 
ACM : CS x sin, ACM= SG. But sin. CMS : SG : : sin. SMN -. 

SN, that is, a : CS x sin. ACM : : s : ^i*- 5 ™' ACM s = SN, and 

therefore ac — bs + CS x si g n - ACM - s = TN + SN = ST = 

CMS — z. Let d— - * s '^'' r- b, and then at -f- ds = 

CMS — #, and CMS = % -f- ac -f- <&> as in the preceding 
method. 

Third Method. 

Let % = the angle SMN, 5 — the series expressing its sine, 
and c == the series expressing its cosine, as before ; but let a 
now denote the sine of CSM, and b its cosine, and let MN 
meet CS in H. Then ac + bs is equal to the sine of the sum 
of the angles SMH, MSC, that is, ac + bs =s sin. MHA = sin. 
ACE, the excentric anomaly. We have therefore CE : ac -f- bs 
: : CS : ™S?!1+M — ST = CMS — %. 

Let a 7 = ^c^~j an( * * = elf- , and then dc + «=r CMS — 2, 

and CMS =5 »4- < fc + « = * + ^-~- + ~.— da " 



2 ' 2.3.4 2.34.5.6 

. r, . es? , est 5 n„ j , , tfc 2 ez z 

+ &C. +*» — — + Ii -^^&C. = £/+*+« - T - - 

JL±Ljl J*! &c Let/ = CMS — d, and then /= Az 

~ 2.3.4 * 2 «3-4-5 ^ ' J 

— B% 2 — C2 3 + Dz* -f Es 5 — &c. putting A, B, C, &c. for the 
coefficients. By reversing this series, or by putting d for a t 
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e for d, and/ for e, in the series in the first method, we find 
rv— H-fJ- R ° d /"J, R ° g f 3 I R ° rf * f 3 - R0< * f* 

I prefer this method to the first or second, and therefore I 
proceed to calculate by it. 

Example I. 

Let us suppose with M. Delambre* that the mean ano- 
maly is 135 , and the excentricity of the orbit 0.25, the mean 
distance from the sun being 1. 

Here CM+CS=i.2 5 , CM-CS= , 75, and §*=g = iIL, 
the log. of which is to be used for any given mean anomaly 
in the orbit. 



•75 

1.2c 



Log. tan. ©^..go' 
Log. tan. 55 ..32.49.84, 



- - Log. 9.7781513 

10.3827757 



10. 1609270 



CSM = 122 ..,52..49.84, 

CMS = .12.. 7.. 10. 16=43630" 16 

206264.8 

CMS = . 2115249 
Log. 9-39794QO 



CMS is found by this 
proportion, 206264^8 : 1 
: : CMS in seconds : its 
length in parts of the ra- 
dius. 



CS = 0.25 
a 



Log. 9 9241783 



£?=. 2099511 Log. 9.3221183 
.2115249 = CMS 

.0015738 =/, Log. 7. 1969495 



CS 

b 



Log* 4-%)/867 
Log, 5.3 i 44251 

Log. 9'3a536i6 

Log. 9-397940O 
Log 9.7347108 



e=.i357 223 Log. 9.1326508. 
As CSM is obtuse, e is negative. 
1 +£=.8642778 Log. 9 9S 66 S34^ 



• Vol. II. page 28. 
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The angle SMN is therefore calculated from the series in 
the following manner : 
First term. 
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R° 

l -{- e 

Numb. 
For 2d term 



Log. 1.7581226 
- Log. 7-1969495 

^■95507^ 
• Log. 9-9366534 



.1043322 

.0000231 



Sum - o°. 1043553 
= SMH. CSM 

ACE =s - 



Log. 9.0184187 
6'..i5."6'7 

= 122°..52..49, 84 



Second term. 



Log. 1.7581226 
Log. 9.3221183 
Log. 4-3938990 

15-4741399 

0.3010300 



R° - 
d 

2, Log. 

15-1731099 
(1+e) 3 - Log. 9.8099602 

Numb. .0000231 Log. 5.3631497 



122..59..5.51 

This differs from M. Delambre's conclusion only in the 
second place of the decimals. 

Example II. 
Supposing, with M. Delambre, that the mean anomaly in 
the same orbit k q6°, required the excentric anomaly. 
We have as before ^X - - - - . Log. 9.7781513 



•75 
.25 



Log. tan. 48 .. 

Log. tan. 33 ■4o / ..4^ / .5i . _ _ . 

CSM = 81 .40 41^,51 
CMS = 14 ..19 ..18. 49 = 51558.49 

206264.8 

CMS = .2499627 



10.0455626 
- 9-8237139 

Log. 4.7123003 
L °g - 5-3H4 2 5 i 
Log- 9-397 8 75a 



*34 

CS = 0.25 
a 
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0.25 



d = . 2473677 

. 2499627 = CMS 

. 0025950 =/ - 



Log. 
Log. 

Log. 



9-39794QO 
99954030 

9-3933430 

7.4141374 



R° 

/ " 

(l + e) 

Num. 
ed term 



Log. 

First term of the series. 

Log. 1.7581226 

L ° g * 7-4141374 

9.1722600 

Log. .0154366 



.1434906 

.0000429 

Sum o°. 1435335 

^SMN. CSM 

ACE 



Log. 9.1568234 



- = 8i c 



8'..36".72 
.40 ..41 .51 



CS 
h 



e = .0361832 



Log. 9.1605669 



L °g- 8.5585069 
1+^=1.0361832 Log. .0154366 



R° 

d 

r 



Second term of the series. 

- Log. 1.7581226 

- Log. 9.3933430 

Log. 4.8282748 



Log. 



!5'97974 4 
.3010300 



(1 + O 8 



15.6787104 

Log. .0463098 



Num. .0000429 Log. 5.6324006 



81 ..49. .18 .23 

This also differs from M. Delambre's conclusion only in 
the second place of decimals. 

Example III. 
Let us suppose ALP to be the orbit of the comet which 
appeared in 1682, and reappeared in 1759, according to the 
prediction of Dr. Halley ; that CE is equal to 18.07575, that 
that CS is equal to 17.49225, and that the mean anomaly is 
1 79°"47'-3 2 "-i7> J* * s required to find the excentric anomaly, 

~ Log. 8.2149815 
12.7415531 



Here 



cm— cs 



.5835 



CM + CS 35.568 

Log. tan. 89°..53'..46".o8 

Lon. tan. 83 ..41 ..34 .5 - - 10.9565346 

CSM = 173 ..35 ..20 58 

CMS= 6 ..12 ..11 58 = 22331.58 Log. 4.3489194 

206264.8 Log. 5-3i44 2 5 5 



CMS 



.1082666 Log. 



9-0344943 



CS = 1749225 
a - 

CE = 18.07575 
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17.49225 



■m 



Log. 1.3428457 
Log. 9.0478932 



d 



10.2907389 
Log. 1.2570963 

.1080544 Log. 9.0336426 
.1082666 = CMS 



11° 
/ 



.0002122 sssf, Log. 6.3^67454 

First term of the series. 

Log. 1.7581226 
- Log. 6.3267454 



(l+e) 

Num. .3171788 



8.0848680 
Log. 8.5835638 



R° 

r 
d 



Log. 9.5013042 

Second term of the series. 

Log. 1.7581226 

Log. 2.6534908 

- Log. 9.0336426 

1 3-445«5 6 o 
Log. .3010300 



(i+£) s 

Num. 



13.1442260 
L °g- 5.7506914 



.0024748 Log. y, 393534 ,6 

.3171788 

.0000386 

.3196922 sum of the positive 
.0000406 



CS 
b 



Log. 1.2428457 
Log. 9-997^757 



11.2401214 
CE=i8. 07575 Log. 1.2570963 

e = . 96166^8 Log. 9.9830251 
As CSM is obtuse, e is negative. 

1-K—0383322 Log. 8.5835638 
Third term of the series. 



R° 

e 

P 
6 

Num. - 



Log. 1.7581226 
Log. 9.9830251 
Log. 8.9802362 

20.7213839 

Log. .7781513 



Log. 



19-9432326 

43342552 



R° 

d % 

P 

2 



— .0000406 Log. 5.6089774 

Fourth term of the series. 

Log. 1.7581226 
Log. 8,0672853 
Log. 8.9802362 

18.8056440 
Log. .3010300 



[i+ef 



Num. .0000386 



18.5046140 
Log. 2.9178190 

L °g- 5-5$6795o 



SMN=o°.3i965i6 = i9'..io".y S 

CSM= 173..35 ..20.58 



1 73-54. -31 -33 = ACE. With the design already 
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expressed, I adopted in this example the same data with Mr, 
Ivory : see Transactions of the Royal Society of Edinburgh, 
Vol. V. page 236. 

The preceding method bears a nearer resemblance to that 
given by Keill, in his Astronomy, than to any other of which 
I know. Adapting his manner of proceeding to the figure here 
used, he putsjy = EM, e = sin. AM,/ cos. AM, and^ = CS. 
The series expressing the sine of AE, is therefore equal to 



e — fy — _i_ JL. jl _£Z_ & c . But the radius, which is 1, is 

JJ 2 » 2.3 ■ 2.3.4 



2 ~ 2.3 "T" 2.3.4 

to the sine of AE as CS or g is to ST or EM, that is to y. 
Consequently y = ge—gfy — ^J- + ^^ + —■ &c. and there- 
fore^ =y+gfy + Sf »^_fg&c. 

By reversing this he obtains a series, which, omitting the 
numbers in the coefficients, converges as the powers of ~-" f 

or 1 + cs x"cos. am * ^ n ' s degree of convergency to the value 

of y in the foregoing examples is as follows. 

In the first example as the powers of .2147372, 
In the second example as the powers of .2553020, 
In the third example as the powers of .1086700. 
In the third method which has been here investigated, the 

series converges to the value of SMH as the powers of 

/ ^ CMS— CS x sin. CSM ,_,, . , r 

"e or 1 + cs X C08 .-CSM • rhls de S ree of convergency in 

the foregoing examples is as follows. 

In the first example as the powers of .0018209, 
In the second example as the powers of .02504,38, 
In the third example as the powers of .0055358. 
The third method, by which the three foregoing examples 
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are calculated, appears to me the most simple and precise in 
theory, and the most expeditious in practice of any which I 
have seen. This I say with the greater degree of freedom, 
as I am so well aware of the similarity between its series and 
that of Keill's, and so perfectly convinced of the advantages 
which it derives from Cassini's approximation, that I consider 
it, with the exception of some deviations, as a combination of 
their methods. 
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